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Abstract 

In the framework of the spectral integral equation, we consider the bb states and 
their radiative transitions. We reconstruct the bb interaction on the basis of data for 
the levels of the bottomonium states with J PC = 1 , ++ , 1 ++ , 2 ++ as well as 
the data for the radiative transitions T(3S) — ► 7Xbj(2P) and T(2S) — ► 7Xbj(lP) with 
J = 0,1,2. We calculate bottomonium levels with the radial quantum numbers n < 6, 
their wave functions and corresponding radiative transitions. The ratios Br[xbj(2P) —> 
7T(25)]/Br[x6j(2P) — > 7T(15)] for J = 0,1,2 are found in the agreement with data. 
We determine the bb component of the photon wave function using the data for the e + e~ 
annihilation, e+e" -> T(9460), T(10023), T(10036), T(10580), T(10865), T(11019), and 
predict partial widths of the two-photon decays rjbo — > 77, Xbo - > 77) Xb2 77 for the 
radial excitation states below BB threshold (n < 3). 



1 Introduction 



In [1], see also [2], the program was suggested for the reconstruction of soft quark-ant iquark 
interaction, on the basis of data on meson levels and meson radiative transitions. Now, we have 
the first results in the realization of this program for the bb, cc systems and light quarkonia qq. 
In this paper we present the results for the bottomonia. 

In [1], the equations for qq systems were written in terms of the spectral integral represen- 
tation. The spectral integration technique is precisely advantageous for composite particles, 
for the content of a composite system is thus strictly controlled and there is no problem with 
the description of high spin states. The equation for the composite qq system in the spectral 
integration technique [1] is a direct generalization of the dispersion N/D equation, when we 
represent the iV-function as a sum of separable vertices. In [1], this equation was conventionally 
called the spectral integral Bethe-Salpeter equation. However, it should be emphasized that 
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in certain important points it differs from the standard Bethe-Salpeter equation [3] written in 
the Feynman technique (the application of the Feynman technique to the calculation of meson 
states may be found, for example, in [4, 5, 6, 7, 8] and references therein). 

The strong QCD is responsible for the formation of composite systems, though, despite a 
remarkable progress in this field, certain features of the soft interaction of color objects are 
rather enigmatic till now. A fascinating feature of the light quark (q = u, d, s) meson spectra is 
the linearity of their trajectories on both the (J, M 2 )- and (n, M 2 )-planes (J is the spin of the qq 
system with the mass M and n is its radial quantum number) [2, 9]. The linearity of trajectories 
has been also seen in the baryon sector [2, 10] — this fact allows one to suggest a specific role of 
diquarks in baryon systems, see the discussion in [2, 10] and references therein. Unfortunately, 
the number of experimentally observed excited states for baryons is much less than that given 
by the standard quark model calculation (e.g., see [11]) that limits our understanding of the 
three-quark interactions. 

In practice, all the observed highly excited meson states (with M > 1500 MeV) are lying 
on linear qq trajectories, and we have no candidates for the qqg hybrids. The effective gluon g 
has the mass m g of the order of m g ~ 700 — 1000 MeV, [12, 13, 14, 15], so the hybrid mesons 
could be spread over the region 1500 — 2000 MeV. Still, the experiment does not point to the 
increase of the meson state density in this region. Also, there are no definite indications to the 
existence of mesons with the exotic quantum numbers inherent in hybrids [16]. 

The particulars of hadron mass spectra discussed above point to the fact that in studying 
mesons in the soft region (or formed at large distances) the most reliable way is to consider 
the quark-antiquark interaction by determining their characteristics from the data. It is the 
guideline in our approach to the quark-antiquark systems. 

The spectral integral equation [1] gives us a unique solution for the quark-antiquark levels 
and their wave functions, provided the interquark interaction is known. Let us emphasize that 
the equations work for both instantaneous interactions (or those of the potential type) and 
the t-channel exchanges with retardation, and even for the energy-dependent interactions: this 
follows from the fact that the equations themselves are the modified dispersion relations for the 
amplitude. For solving the inverse problem, that is, for reconstructing the interaction, it is not 
enough to know the meson masses — one should know the wave functions of quark-antiquark 
systems. Such an information is contained in the hadronic form factors and partial widths of 
radiative decays. Therefore, in the present approach, we consider simultaneously the meson 
spectra in terms of the spectral integral equations and meson radiative transitions in terms of 
the dispersion relations over meson masses — in this way, all the calculations are carried out 
within compatible methods. 

The method of calculation of the radiative transition amplitudes in terms of the double 
dispersive integrals was developed in a number of papers [17, 18, 19]. An important point was 
the representation of the transition amplitude in the form convenient for simultaneous fitting 
to the spectral integral equation — it was done in [20, 21]. 

A significant information on the quark-antiquark meson wave functions is hidden in the 
two-photon meson decays: meson — > 77. For the calculation of such processes, one needs to 
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know the quark wave function of the photon; the method of reconstruction of the 7 — > qq and 
7 — > cc vertices was developed in [22, 23]. 

The paper is organized as follows. 

In Section 2, we present the technique and basic formulas which were used in fitting to 
quark-ant iquark states. Here, we briefly recall the spectral integral equation and give the 
formulas for the radiative transition amplitudes of quark systems calculated within the double 
spectral integrals. Keeping in mind the application of these formulas to the other quark systems, 
we do not specify the flavor of the considered quark: (QQ)i n — > 7 + (QQ)out , e + e~ — > V(QQ) 
and QQ-meson — > 77. 

In Section 3, we consider the bottomonium systems: the choice of the bb- mesons as a primary 
object for the study is motivated by the large 6-quark mass, so the nonrelativistic approximation 
is expected to work well, and we can reliably compare our results with those obtained within 
nonrelativistic approaches [24, 25, 26, 27]. Different variants of the interaction (instantaneous 
and retarded) are discussed. The results of fitting to masses and radiative transitions (for both 
observed and predicted states) are presented and discussed. In the bb sector, we reconstruct 
the interaction, on the basis of data on the levels of bottomonia with J PC = 0~ + , 1 , ++ , 
1++ 2++ and radiative transitions T(3S) -> 7X 6 j(2P) and T(2S) -»• 7X 6 j(lP) at J = 0, 1,2. 
We calculate the bottomonium levels with radial quantum numbers n < 6, as well as their 
wave functions and corresponding radiative transitions. We determine the bb component of 
the photon wave function on the basis of data for the e + e~ annihilation reactions e + e~ — > 
T(9460),T(10023),T(10036), T(10580),T(10865),T(11019), and predict partial widths of the 
two-photon decays i] b0 — > 77, Xbo ~^ 77; Xb2 ~^ 77 for the excited states with n < 3. 

Brief summary is given in Conclusion. 

2 Spectral integral equation and radiative transition am- 
plitudes 

Here we present the formulas used in fitting to the bb systems. Of course, they may be also 
applied for the charmonia, as well as for light quark states qq with I — 1, or one-flavor states 
with 1 = (pure ss or nn = (uu + dd) / y/2 systems. 

2.1 Spectral integral equation 

The wave function of the quark-ant iquark meson with the mass M is characterized by the total 
momentum J, quark-antiquark spin S (in the flavor nonet with fixed J p , the values S — 0, 1 
determine C parity) and radial number n. We denote the wave function as ^\ v } ) ^ 1 ... flJ (k±), with 
k± being relative quark momentum and the indices /ii,... , /ij related to the total momentum. 
For the heavy-quark QQ system, the spectral integral equation reads [1]: 

( S - M2 )CN= (1) 
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Am 2 

Here, the quarks are mass-on-shell, k\ = kf —k\ — k' 2 2 = m 2 . The phase space factor in the 
intermediate state is determined as follows: 



d<Z> 2 (P';k[,k' 2 ) = - 



1 d 3 k[ d 3 k' 2 



2 (2?r) 3 2fci (2tt) 3 2^ 



-^^(p'-K-k',) . 



The following notations are used: 

k ± = l - (h - k 2 ) , P = k 1 + k 2 , k' ± = l - (k[ - k' 2 ) , P' = k[ + k' 2 , 



P=s, P' 2 = s' , <?i = 9iXV - 



P P 



V- v „'-L 



9au ~ 9fj.f 



pi pi 



(2) 



(3) 



so one can write k^ = k u g^ and k'^ = k' v g'^ . In the center-of-mass system, the integration 
can be re-written as 



Am 2 



(4) 



where k' is the momentum of one of the quarks. 

For the fermion-antifermion system with definite J, S and L (angular momentum), we 
introduce the moment operators Q^'.b^] (k±) defined as follows [28]: 



«i J (^) = ^ 1 ..., J (^) 



'2J+1 



aj 



(2J+ I) J 



^(J + l) aj ' 



QK:i" /) (*±) = 7«^ 



'J + l 



a/il— Mj 



J 



0K J) (*i) = 7.^Jr. 



where aj = (2J - 1)!!/J! and 



k 2 



J! 



jU-L jU-L jU-L 1.-L JU-L _ 



2J 



1 {pti^ti^ ■ ■ ■ ^ + g^k^k^ . . . k^ + . . .) + 

_*i ( 

- 1) 2J-3) v 



rt-L fl-L i.- 1 - I.- 1 - i.- 1 - 1 

_ 1)(2J — 3) rW 2y ' 13 w ^ « ' ' ' w 
dm^dnsnsk^k^ ■ ■ ■ k^ + ...) + ... 



(5) 



(6) 



(7) 
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2 y q-L x^- 1 ) r£V 

2J _ ]_ y MiMj ;ii...Aii_i/ii+i...;i :7 -_i/i J - + i.../ijaV -U 



The operators are normalized as: 
dVL 



J ^MQ { °; J .il(rn + h)Q^ J :P L {m - k 2 )\ = -2sk 2J (-1) J 0%;;% , 
j ^Sp[Q&5?> + k)Q^ J .i]{rn - k 2 )\ = -2sk^ \-iy '0%;;% , 
/ ^Sp[Q^y)(m + kM^fim - k 2 )\ = C iJ 2 jff ~ 2s)k^\-iyOZl 
J ^MQ^/Krn + kMd^irn - k 2 )\ = - 2 S )k 2 ^)(-l)^ 



V\...VJ ' 

/ ^Sp[Q^.-^)(m + kMd^irn - fc)] = -8 ^ ^ k^\-l) J O^f . (8) 

Here, O^'/'/f™ is the projection operator to a state with the momentum J and s = Am 2 + Ak 2 . 
The projection operators have the following properties: 

x£lv L (k ± )0%;;% = ^(fcj , 0£;;;£ = . (9) 

Let us underline that (— 1) J 0^ i 1 '"/f J describes the spin structure of the propagator of a particle 
with the total spin J (see [28] for more detail). 

In terms of these operators, the wave functions read: 

=Q^(*±)tf' £=J ' J) (*i) , S = 0,landJ = L, (10) 

where the functions ^^' L,J \k 2 1 ) depend on k\ only (recall that in the center-of-mass system 
k\ = -k 2 ). 

The wave functions with L = J are normalized as follows: 

1 -/(^ 2s i k i 2J ^ J ' J, (^i 2 ' ' u > 

while for L = J ± 1 the normalization reads: 



+ 16 V 2J + + i 1) k2(J+lV " 5 '^W^^) + lC' J - M ^l)lX 2 *-^) k2(J ~ 1) 
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We re-write this normalization condition as follows: 

Wj+M+i + + = 1, (13) 

where Wj+ij+i, Wj+ij-i and Wj_i jt /_i are determined by the wave function convolutions 

^{s,j+i,j)^(s,j+i,j)^ ^{S,j+i,j)^{s,j-i,j)} and ^(s,j-i,j)^(s,j-i,j)^ corr e S pondingly. 

The interaction block can be expanded in a series with respect to a full set of the t-channel 
operators Oj\ 

6 1 = I, ip, icr^, ^75, 75 , 
V (s, s', (k ± k' ± )) = Vi (s, s', (k ± k' ± )) O i ®O i , (14) 

By solving the equation, the t-channel operators should be transformed in the two s-channel 
ones; this procedure is clarified in Appendix I. 

The equation (1) is written in the momentum representation, and we solve it in the momen- 
tum representation too. The equation (1) allows one to use as an interaction the instantaneous 
approximation, or take into account the retardation effects. In the instantaneous approximation 
one has 

V (s, s', (k ± k' ± )) — V(t ± ), t ± = (k 1± - k' 1± ),(-k 2± + k' 2± ), . (15) 

The retardation effects are taken into account, when the momentum transfer squared t in the 
interaction block depends on the time components of the quark momentum (for more detail see 
Section 2.5 of [1] and the discussion in: [29, 30, 31, 32]): 

V (s, s', (k ± k' ± )) — V(t), t = (k 1 - k/M-k 2 + k' 2 ), . (16) 

Fitting to quark- ant iquark states, we use the interaction blocks with the following t-dependence: 

47T 



I-l 

Io 



8"7T/Z 



h = 8tt 



(JJL* - t) 2 ' 



> 2 -t) 3 (fi 2 -t)\ 

I 2 = 967T/1 



( 2/i 2 1 



T Qfi I 16 ^ 12/i 2 , 1 \ 

= X> + St) N+1 - n (n - VtT ■ (18) 

\^ l ) n=0 

Traditionally, the interaction of heavy quarks in the instantaneous approximation is represented 
in terms of the potential V(r). The form of the potential can be obtained with the Fourier 



or in general case 
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transform of (17) in the center-of-mass system. Thus, we have 

t± = ~(k - k>) 2 = -q 2 , 
4rW) = /^e^/A), (19) 

that gives 



r (coor) 



N 



(r, //) = e-» r . (20) 



Working with the instantaneous interaction, we consider the following types of V(r): 



^(r) = a + b r + c e'^ r + d , (21) 

r 

where the constant and linear (confinement) terms read: 

j(coor) i p,\ 
a ~~ > a I V, constant ~^ U J , 

br - & Z} ^ (r, Wi»«,r ->0) . (22) 

The limits ^constant f^iinear mean that in the fitting procedure the parameters fi cons tant and 
^linear are chosen to be small enough, of the order of 1-10 MeV. It was checked that the solution 
for the states with n < 6 is stable, when fi CO nstant and Hu ne ar change in this interval. 



2.1.1 Wave functions 



Let us present in the explicit form the wave functions of the studying states: + , 1 , ++ , 
1 ++ , 2 ++ , l +_ . The wave functions of these states, ^[f)^ M2 in terms of the operators (5) 
and invariant functions ^^^'^(k^) read: 



o-+ 

1— 

0++ 

1++ 

2++ 



(n) 

*Si?i=^'°' 1) (* 2 ) + ^ 



A1) (* 2 ) , 



(n) 
(n)/x 



mfy#' ll0) (fc 2 ) , 
3 



IE 



€^\k 2 ) , 



+ 



(1,2) 

(n) /HI ^42 

5 

,(o,i) 



kixil^ + kn.^^ rkg^ lfi2 



^n lX2 \k 2 ) + 



^• 3 ' 2) (A; 2 ), 



*g) 1 i = V3<75M?' 1>1) (* 2 ) 



(23) 



Here, the short notations are used k 1 - = k = (hi — k 2 )/2 {(kP) = at quark masses equal to 
each other) and £ 7 p^ = e UlU2U:ifl / -f Ul P U2 k u . i as well as the equality k = ml (recall again that in 
the spectral integral technique the constituents are mass-on-shell). 
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The 1 and 2 ++ states are defined by two wave functions with different L; correspondingly, 
there are two wave functions for these two levels. 

The 1 state is characterized by two angular momenta, L = and L = 2. The states in 
(23) are not pure with respect to L, for the wave function r j^ ( - 1 ' ' 1 \k 2 ) is a mixture of S- and 
P-waves. The pure S state is given by the following operator (e.g., see [17, 21]): 



= it 



k.. 



2m + y/s 



(24) 



r ~ 

1 M — 



In the nonrelativistic limit, the operators and 7^7 coincide: 
pure P state of the system 2 ++ is defined by the following operator [21]: 



nonrel 



Likewise, the 



T, 



Mi ^2 



1*2*- £H 



:(kr)g 



M1M2 



(25) 



Note that for heavy quarks the use of pure operators (24), (25) results in the insignificant 
modification of the wave function representation, for the 77} operator leads to a small ad- 
mixture of the P-wave and, vice versa, the ^= kjt — \k 2 ^ operator results in a state with 
dominant P-wave, with a small admixture of the S'-wave. Similarly, for the state 2 ++ the 



operator 



^mi7 M2 + ^M2 7 w 



5 

V2 



IK, 

+ *W7, 



leads to the dominant P-wave, while the operator 
gives us the dominant F-wave. 



in 



Let us note that dealing with the operators of pure states, like (24) and (25), would not 
facilitate the fitting procedure, for real states are the mixture of different waves (S, D for 1~" 
and P, F for 2++). 



2.2 Radiative transitions (QQ)i n — > r f(QQ) ut 

Here, we list the formulas for the amplitudes and partial widths of the radiative transitions 
{QQ)in ~^ l{QQ)out used in the fit. The technique for the calculation of the transition ampli- 
tudes (QQ)in — > l{QQ)out was developed in [18, 19, 20, 21]. In [33], the transition form factors 
were transformed to the form convenient for the fitting procedure, and we present them in this 
form below. 



2.2.1 Transitions of the vector (1 ) and pseudoscalar (0 + ) mesons 

Here, we present the formulas for the radiative transitions of the vector (V) and pseudoscalar 
(P) mesons. The partial widths for the decays V — > 7P and P — > read: 



Mr i, m-Mi? 2 

MV \ (Ml-M v f 2 
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where a = e 2 /4n = 1/137. The transition form factor is expressed in terms of the QQ wave 
functions ^°>°\k 2 ) and ^ 1,0>1) (A; 2 ), see (23). However, in line with [19, 20, 21], we change the 
notations as follows: 



In terms of Vv( s ) an d ipp(s), the form factors are written as: 

F V(0) ^ P = Z v ^p— \ — ip v{0) (s)ip P (s)ki — 

4:71 J 7T x /s — \/s — 4m z 



V(2)^ 7 P 



z 



4m 2 



y-> 7 p 



32V/2 7T 



4m 2 



7T 



^V(2)(s)^p(s) X 



(27) 

(28) 
(29) 



x 



y/s + Vs — 4m 2 



(2m 2 + s) In ^ - y 7' ^= - 3^/s(s - 4m 2 ) 



y/s — y/s — Am 2 



where Zy^~ jP is the charge factor: 

Zy-> 7 p = 2eQ, 
The total transition form factor reads: 



e c - 3 , e b - 3 . 



Fy^ y p — Fv( )^ 7 P + i 7 V(2)-» 7 P • 



(30) 



(31) 



Normalization conditions are determined by Eq. (12), they can be written, after the integration 
over the angle variables, as the integrals over s: 



r d 3 k r 

J (27T 3 )k J 

Am 2 



ds / s — 4m 2 



16tt 2 



Normalization conditions for vector state determined by Eq. (12) read: 



W 00 [V}+W 02 [V}+W 22 [V} : 



W 00 [ V] 
W 02 [V] 
W 22 [V] 



1 

3 



4m 2 



ds 
16t? 



^ (0) (s)4(s + 2m 2 ) 



4m 2 



(32) 

(33) 
(34) 



/ ^v(o)(s)VV(2)(s) (s - 4m 2 )' 



' s — 4m 2 



4m 2 

oo 



ds 



3 / iS 2 " ^ (s) 



Sm 2 + s)(s - 4m 2 ) 2 /s-4m 2 



4m 2 
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Normalization condition for ipp(s) in the s-integral representation is: 

1 



ijjp(s) 2s 



' s — 4m 2 



4m 2 



16tt 2 



(35) 
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2.2.2 Transitions of the vector (1 ) and scalar (0 ++ ) mesons 

The partial widths of the vector (V) and scalar (S) meson decays, V — > 7S' and S — > 7 V, read: 

1 Afy — Af| 2 

Myry^ 7 5 = -a ^ |-fV-»7S| ' 

1 Ml — M?r . |2 

MsTs^v = -a g ^ |^ 7 5| 2 - (36) 

The transition form factor Fy^s is expressed in terms of the wave functions ip£' 0,1 \k 2 ), 
*Pn' 2 ' 1] ( k2 ) and Vi 1 ' 1 ' ^^ 2 )- Changing the notation in line with, [19, 20, 21] 

^ l > \k 2 )=Ms), (37) 



we have: 



2 w 

i*V(o)- 7 s = Zv ^ s ^ j -^i'v(o)( s )i's(s)I v ^ 1 s(s), 



4 m 2 

y/2,fn 2 f ds 

F V (2)^s = Z v ^ys 16?r y -^(2)(s)^s(s)(-s + 4m 2 )/i/_ 7S (s), 

4m 2 



y/s + Vs — 4m 2 



Iv^s(s) = Js(s - 4m 2 ) - 2m 2 In v I ' v "r= ^= ■ (38) 

V s ~ V s — 4m z 

The total form factor is equal to 

Fv-*-iS = F V (o)^-ys + ^V(2)-t- 7 s • (39) 

The charge factor of the transition V — > 76" for heavy quarks coincides with that for V — > 7P: 
-^v->7S — -^v-^7P? see Eq. (30). Normalization condition for ips( s ) reads: 



1 = / i^ ^( s ) 2m2 ( s - 4m2 )V iL T^- (40) 

4m 2 

Working with the transition amplitude V — > jS, one faces the ambiguity in writing the spin 
operator that is due to the existence of the nilpotent spin operator component — this problem 
is discussed in detail in [34, 35]. 

2.2.3 Transitions of the tensor (2 ++ ) and vector (1 ) mesons 

Three independent spin operators determine the transition amplitude T — > 7V and, corre- 
spondingly, we have three form factors [33]. In terms of the form factors F^, v (i = 1,2,3), 
the partial widths for the decays T — > 7V and V — > 7T read: 

m V - a m T- m V L± (p(l) N2, l/ p (2) N2, l/p(3) x 2 

ZAj Tflrp 

m V - a m V- m T [ I (p (X) x2,_±/ p (2) x2,_±/ p (3) x 2 l N 
m y ly^ 7T - — —2 Z n (-^ T ^ 7 y J +^22l^T-» 7 vJ + Z 33 i^T^v) ■ l 4i J 



v 
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Here, 



3M| + MMj,My + 3M V 



-22 



UM^My 
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M| + lOM^My + My 

3.U/ • 3I.U7.U, 2 • 3.U, 1 ' 



Z 33 — ~ 



9 (M 2 + M 2 ^ 2 



2 Mt + 10M 2 M 2 + Mt ' 



(42) 



L v 



The quark-ant iquark 2 ++ state is determined by two components of the wave function with 
the dominant P- and F-waves: 



With these notations, the form factors -^r(L)^ 7 y(L') for i = 1, 2, 3 read: 



(43) 



Here, 



where 



^T(LH 7 V(L') - ^T(L)-> 7 y(L') 



^T(1H 7 V(0) 
c(2) / 

°T(i)^ 7 y(o) ^ 6 



,(i) 

'T(l 
7 (2) 

'T(1)-»7V(2) 
'T(3)-»7V(0) 



St(i)->jV(2) 1 6 



;(2) 
'T(3 

7 (1) 



^T(3)^ 7 y(0) 1 * 



^T(3)^ 7 y(2)l S 

C(2) / 
°T(3)^ 7 F(2) 



4Vw 

41rv(«) 



4m 2 



16^ 



^T(L)- +7 y(L')( s )V ; r(s)^(s) . 



-^(8m 2 + 3 S )4V( s )- 
2 c(3) ^ 



4i 7 v( s ) ) 



40 

7 (3) 



(16m 2 -3s)(4m 2 -s)4 J 
V2 



(i) 



3 ^/^ .2 „\2r(l) 



20 



(3) 



i4m'- s yi^ v (s), 

>) = -^(6m 2 + S )// 2) 
(4m 2 -s) 2 (8m 2 + 5)4^(5) , 



3 ^T(3)-> 7 y(o) 
3 

~80 

2 o(3) 

3^T(3)^ 7 y(2)> 



(s) = -^(16m 2 - 3s)(4m 2 - S )41 yV (s 



72 



2m In 



y^i + Vs — 4m 2 
^/s — Vs — 4m 2 



s(s — 4m 2 ), 



m 2 (m 2 + s) In ^ + V § ^ - ±- J S (s - 4m 2 ) (s + 26m 2 ). 
V^s — v s — 4m 2 12 v 



(44) 



(45) 



(46) 
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The total form factor is a sum over four terms: 



Normalization condition for the tensor meson wave function is written in the s-integral repre- 
sentation as follows: 

1 = W U [T] + W 13 [T] + W 33 [T}, (48) 

Wn[T] = I / ^(8m 2 + 3 S )( S -4m 2 )^i^, (49) 

4m 2 
oo 



TT , 1 f ds / x , / n v^3 . r,., s — Am 2 



4m 2 

oo 



^ 33 [T] = ^/ T ^^(3)(^)^(6- 2 + S )( S -4m 2 ) 3 ./ S ^ 

4m 2 



2.2.4 Transitions of the pseudovector (1 ++ ) and vector (1 ) mesons 

Here we present the formulas for the radiative transitions of the vector (V) and pseudovector 
(A) mesons. The partial widths for the decays A — > 7V and V — > are determined by the 
form factor Fa-^v as follows [33]: 

a m\-m\ , 2 

«U1 A^yV — T7) ~2 2 AV ^A^yV > 

^r^ 7 A = ^ 2 a 4aF 2 a ^ v , (so) 

where 

± Ml + 6M 2 M 2 + My ± M\ + 6M 2 M 2 + M v 

Zav ~ 2Mj ' ZyA ~ 2M| • (5ij 



Changing notation 
we write: 



€ Xl \k 2 )=Ms) , (52) 



/ds 
Y^S A ^v(L){s)lpA{s)lpV(L){s) , 

Am 2 

S A ^v(o)(s) = -]J^Ia%v( s ), 

*W<2)(*) = ^(4m 2 - S )/iV( s )' ( 53 ) 
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where 



/£> (s) = f 2m 2 In ^ + V / - Js(s - 4m 2 ) V (54) 



The normalization condition reads: 



y^i — Vs — 4m 2 



1 = i / *(* " 4- 2 ) \ 8 —^- (55) 



2 7 16tt 2 

4m. 2 

2.3 Radiative transitions e + e~ — > and QQ-meson —> 77 

For the consideration of the radiative transition processes e + e~ — > V(QQ) and 77 — * QQ- 
meson, it is convenient to introduce the quark-ant iquark components of the photon wave 
function [19, 22, 23]. The quark components of photon is tightly related to the determination 
of the quark wave functions of vector mesons. 

The introduction of the quark-ant iquark photon wave function may be illustrated by the 
two-photon meson decay. Dealing with the time-ordered processes, that is necessary in the 
dispersion relation or light-cone variable approaches, the QQ-meson — > 77 decay should be 
treated as a two-step reaction: the emission of photon by the quark (Fig. la) or antiquark 
(Fig. lb) and a subsequent annihilation QQ — > 7. The triangle diagram cuttings related to 
these two subprocesses are shown in Fig. lc, thus leading to the representation of the triangle 
diagram in terms of the double dispersion integral. In the diagram lc, on the left from the first 
cutting, there is the transition vertex of quarkonium — > QQ: we denote this vertex as Gqq(s). 
This vertex determines the wave function of the initial QQ-meson: 

^ = *»(»)■ < 56 > 

In the previous subsection, this wave function was denoted for different J PC as Vv(s), i>p( s ), 
ijjs(s), and so on. 

Likewise, the right-hand cut in Fig. lc describes the transition QQ — > 7 and provides us 
with the factor 

^7 e Q » (57) 

where s' is the invariant energy square of quarks in the final state and cq is the charge of the 
Q-quark. But when we deal with the transition QQ — > 7, the interaction of quarks should be 
necessarily taken into consideration. 

The quarks may interact both in initial (Fig. Id) and final (Fig. le) states. In fact, 
the interaction of quarks in the initial states has been accounted for in (56), because the 
vertex functions Gqq (or wave functions tpp, ips, an d so on) are the solutions of the spectral 
integral equation — this equation is shown diagrammatically in Fig. If. As concerns the 
quark interaction in the final state, it should be specially taken into account, in addition to the 
pointlike interaction (57). The diagram shown in Fig. lg stands for the quark interaction in 
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the transition QQ — > 7, and we approximate it with the sum of pole terms of the vector meson 
(T's or ^'s in the cases of bb or cc systems), see Fig. lh. Accordingly, the factor related to the 
right-hand cut of Fig. lc is written as follows: 



e Q , (58) 



s' 

where the vertex function G^qq^s') at s' ~ 4mg is the superposition of vertices of the V{n)- 
mesons (see Fig. lh): 

G^ q q{s) ~ C n G V {n){s) , s ~ Am% ■ ( 59 ) 

n 

Here, n is the radial quantum number of V-meson and C n 's are the coefficients which should 
be determined in the fit. 

At large s, the vertex QQ — > 7 is a pointlike one: 

G^qq(s) ~ 1 at s > s . (60) 

The parameter sq can be determined from the data on the e + e~ -annihilation into hadrons: it 
defines the energy range where the ratio R(s) = a(e + e~ — > hadrons) / a (e + e~ — > reaches 
a constant-behavior regime above the threshold of the production of heavy mesons. The data 
[36] give us s ~ (10 - 15) GeV 2 for the cc component and s ~ (100 - 150) GeV 2 for the bb 
one. 

Therefore, to describe the transition QQ — > 7 we may introduce the characteristics which, 
similarly to (56), can be called the QQ- com P onen t, of the photon wave function: 

= W) ■ («D 

Here g is the photon four-momentum. Let us emphasize that such a wave function is determined 
at s > Arrig. 

There are more reactions which are promptly determined by the photon wave function. 
These are the transitions e + e" — > T and e + e~ — > ip, see Fig. 2: here the loop diagram is 
defined by the convolution of the vector meson wave function and G^qq. 



The transition 7 — > QQ is determined by two spin structures, 7 Q and | k a k — ^k 7^ 
Eq.(23)) and, correspondingly, by two vertices. 



sec 



^l QQ {s) , 7*4? G £W«) (62) 
It means that we take into account the normal quark-photon interaction, j a , as well as the 
contribution of the anomalous magnetic moment. 

For the vertex function of the transition 7 — ■> QQ we use the following fitting formula: 

<2««M = E Wn.«,W + 1 + exp( _^ (s _ So)) • (63) 

6 

G ^Qq( S ) = 2 C nD Gv(nD)(s) , 
n=l 
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where G v{nS) (s) = ^ 101) (s)(s-M^ (riS) ) and G v{nD) (s) = ^ 2l \s)(s-M 2 {nD) ). The parameters 
C n s, C n c, f3 7 and s are determined in the fit. 

2.3.1 Decay V -> e+e 

The transition amplitudes V — > e + e~ were calculated in [23, 37]. The partial width for the 

decay V — > e + e~ reads: 



r(V -> e+e - ) = 



7Ttt 2 



My — 4u 2 /8 „ 4 9 \ , ,o 

1 ' e (g^+gM^jlF^^-l 2 , (64) 



where fi e is the electron mass, My is the measured vector quarkonium mass, and the form factor 
Fy^e+e- is determined by the process of Fig. 2. The quark loops in Fig. 2 are different for the 
S- and D-wave QQ states, so we have two transition amplitudes i*V(n,L)-t-e+e- , with L = and 
L = 2: 

^V-^e+e- = i 7 V(0)- + e+e- + -^V(2)- + e+e- , (65) 

where 



F v ^ e+e - = Z^y[W c J ^v m (s)G\% Q (s) (Jm 2 + A -s) ^-J^L , (66) 

4m 2 

„ „ v ^+*.- ATT f ds , , ,^(n) , As — Am 2 ) 2 l s — Am 2 

4m 2 

At L = 0,2, the wave functions ?/V(l)( s ) are normalized according to (33); here the charge 
factors are Z^ e+e ~ = -1/3 and Z^ e+e ' = 2/3. 

2.3.2 Decay P -> 77 

The two-photon decays of pseudoscalar mesons (L = 0) were studied in [19, 22]. Partial 
width for the decay P — > 77 reads: 

r(P^77) = ^ 2 M||F P _ >77 | 2 . (67) 

The transition amplitude is determined by the processes of Figs, la, 6, it reads: 

ds , , w , . , a/s + Vs — 4m 2 



F P ^ = ZJ^m / ^H.)*^ W >« ' <68) 

4m 2 V 

Recall that ^ 7 _ + qq(s) = G 1 _qq{s)/ s. To underline the existence of decays with different radial 
excited states, we introduce the index n in (68). Normalization of ipp(s) is given by (35), and 

Z ll = I ZS = I (69) 
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2.3.3 Decay S -> 77 



The two-photon decay of the ++ qq mesons (L = 1) was considered in [19, 20]. The partial 
width of the decay S(n) — > 77 (n is radial quantum number) reads: 



T(S - 77) = 



7TCT 



Mc 



(70) 



x 



with the quark transition amplitude (Figs. la,b) equal to: 

1 — f ds 

Am 2 

x Js(s - Am 2 ) - 2m 2 In . 

\ Vs — v s — 4m 2 / 

Normalization of ips{n)( s ) is given by (40). 



(71) 



2.3.4 Decay T -> 77 

The two-photon tensor meson decay amplitude was calculated in [19, 21]. The partial width for 
the decay process T — > 77 is defined by two transition amplitudes with the helicities H = 0, 2: 



r(T - 77) = 7 



4 7rcr 



5 Mn 



7 (0) 
T-»77 



+ 



7 (2) 
T^77 



(72) 



Taking into account the P- and F-wave quark-ant iquark component, we write the form factors 

(73) 



42ry as 



r T^~ r/ - - r T(l)^77 ^ - r T(3)^77 



where 



/.•(" ) _ 777 



4m 2 



ds 
16T 2 



^T(L)(s)^qq(s)4Tl)^77( S ) , 



(74) 



with the following spin factors 

4 



? (0) 

5 T(3)-77 



» = 



2 x /2s (s - 4m 2 ) 



(72m 4 + 8m 2 s + s 2 ) + 



(75) 



m 2 (8m 4 + 4m 2 s + s 2 ) In 



s + y/s (s — 4m 2 ) 
s (s — 4m 2 ) 
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and 



%^ 77 («) = 3^3 (5m + s) - (2m + s) In ■ 



s — y s (s — 4m 2 ) 

4| H77 ( S ) = -V-J_ L ( 30 m 4 - 4m 2 S + S 2 ) - (76) 

2v/2 2 , 4 2 2 n s + v/s (s-4m 2 ) 
-m (12m -2m s + s lln 



5 s — ^ s (s — 4m 2 ) 

Normalization of V ; t(i)(s) and ipT(3)(s) is determined by (48). 



3 Bottomonium states found from spectral integral equa- 
tion and radiative transitions 

The quarkonium wave functions are fitted in the following form: 

^ J \k 2 ) = e-? k2 j2 Ci (S,L, J;n)^-\ (77) 

i=l 

where we re-denoted k 2 = k 2 ; recall that s = 4m 2 + 4k 2 . The fitting parameter (3 is of the 
order of 0.5 — 1.5 GeV~ 2 and may be different for different flavor sectors. We put here /3 = 1.2 
GeV- 2 

The data in the bb sector can be described by two types of the t-channel exchanges, scalar 
and vector ones: 1(8)1, 7 M ®7 M . The addition of the pseudoscalar exchanges like 75 (8)75 does 
not improve the fit. 

Here, we present three variants of fit: with instantaneous forces (solution I(bb)), retarded 
interactions (solution R(bb)) and that with the universal "confinement potential" - instanta- 
neous interaction with nearly the same parameters a and b in (21) for all the quark sectors: bb, 
cc and qq (solution U(bb)) . 

For the bb sector, the parameters for scalar and vector exchange interactions (I (8 I and 
In <8 J ^ see Section 2.2) are as follows (all values are in GeV): 



Interaction 


Wave 


a 


b 


c 


He 


d 






I(bb) 


-0.151 


0.160 


0.506 


0.201 


-0.250 


0.201 


(I® I) 


U(bb) 


0.911 


0.150 


-0.377 


0.401 


-0.201 


0.401 




R(bb) 


-0.680 


0.130 


1.322 


0.201 


-0.228 


0.401 




I(bb) 


-0.812 


0.000 


0.867 


0.401 


0.300 


0.001 


{in ® In) 


U(bb) 


1.178 


-0.150 


-1.356 


0.201 


0.500 


0.001 




R(bb) 


-1.620 


-0.005 


1.821 


0.201 


0.311 


0.001 
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3.1 Masses of the bb states 



The fitting procedure prefers for the constituent 6-quark the mass value m b = 4.5 GeV. This 
value looks quite reasonable if we take into account that mass difference of the constituent and 
QCD quarks is of the order of 200 — 350 MeV and the QCD estimates [38] give the constraint 
4.0 < m h {QCD) < 4.5 GeV. 

The masses of bb states for n = 1,2,3,4,5,6 (experimental values and those obtained in 
the fit) are given below, in (79) - (84). The bold numbers stand for the masses, which are 
included in the fitting procedure. In parentheses we show the dominant wave for bb state (S or 
D for 1 and P or F for 2 ++ ). The right-hand side columns show the mean square radii of 
bottomonia (in GeV" 2 for solution I(bb). 

We have the following masses (in GeV) for 1 states: 



for + states: 



for ++ states: 



otate 


Data 


I{bb) 




U(bb) 


R(bb) 




D2 

Rj 


I [i-Oj 


y.4ou 


9.392 (S) 


9.382 (S) 


9.448 


(S) 


n o a o 
U.o4z 


I [zb j 


lU.UZo 


10.029 


(S) 


10.027 (S) 


10.023 (S) 


1 con 


T(LD 


) 10.150 


10.159 


(D) 


10.158 (D) 


10.105 (D) 


0.342 


T(3S) 


10.355 


10.368 


(S) 


10.365 (S) 


10.362 (S) 


3.794 


T(2D 


) 10.450 


10.439 


(D) 


10.436 (D) 


10.156 (D) 


1.632 


T(AS) 


10.580 


10.615 


(S) 


10.634 (S) 


10.62 


8(S) 


6.504 


T(3D 


) 10.700 


10.661 


(D) 


10.677 (D) 


10.430 (D) 


3.794 


T(5S) 


10.865 


10.819 


(S) 


10.872 [s) 


10.851 (S) 


9.793 


T(4L> 


) 10.950 


10.852 


(D) 


10.898 (D) 


10.669 (D) 


6.504 


T(6S) 


11.020 


11.019 


(S) 


11.084 (S) 


11.040 (S) 


11.990 


T(5D 




11.023 


(D) 


11.109 (D) 


10.892 (D) 


9.793 


T(QD 




11.214 


(D) 


11.303 (D) 


11.085 (D) 


11.990 , 




State 


Data I(bb) 


U{bb) R(bb) 


Rj 






m{is) 


9.300 9.334 


9.322 9.393 


0.922 








10.006 


10.011 10.007 


2.782 






Vb(3S) 


10.344 


10.355 10.352 


5.781 








10.557 


10.626 10.621 


18.839 






Vb(5S) 


10.636 


10.864 10.846 


13.699 






Vb(6S) 


10.837 


11.079 11.037 


11.668 , 




State 


Data 


I(bb) 


U{bb) 


R(bb) 


Rj 






Xbo(lP) 


9.859 


9.852 


9.862 


9.851 


0.847 






Xbo(2P) 


10.232 


10.241 10.236 


10.227 


2.632 






X & o(3P) 




10.509 10.517 


10.510 


5.161 






X & o(4P) 




10.726 10.759 


10.745 


8.053 






Xbo(5P) 




10.88 


4 10.983 


10.952 


12.437 




Xbo(QP) 




10.947 11.185 


11.084 


19.969 , 



(79) 



(80) 



(81) 
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for 1 ++ states: 





Data 


T(hh\ 

1 (00) 




1 h I UU J 




y M f IP) 
A MA- 1 -- 1 / 


9.892 


9.884 


9.895 


9.890 


0.915 


X6i(2P) 


10.255 


10.257 


10.252 


10.249 


2.777 


Xm(3P) 




10.516 


10.528 


10.526 


5.814 


Xm(4P) 




10.697 


10.767 


10.762 


18.944 


Xm(5P) 




10.759 


10.989 


10.970 


13.544 


Xm(6P) 




10.920 


11.191 


11.199 


11.702 , 



(82) 



for 2 ++ states: 



State 


Data 


I(bb) 




U{bb) 


R(bb) 




Rj 


Xm{lP) 


9.912 


9.909 (P) 


9.911 (P) 


9.923 (P) 


0.956 




10.268 


10.270 


(P) 


10.262 (P) 


10.268 


(P) 


2.782 


Xb2(lF) 




10.365 


(P) 


10.347 (P) 


10.346 


(P) 


0.956 


rfP) 




10.528 


(P) 


10.535 (P) 


10.538 


(P) 


5.361 


Xb2(2F) 




10.594 


(P) 


10.592 (P) 


10.592 


(P) 


2.782 


Xb2^P) 




10.738 


(P) 


10.773 (P) 


10.768 


(P) 


8.573 


rfP) 




10.788 


(P) 


10.813 (P) 


10.807 


(P) 


5.361 


rfP) 




10.846 


(P) 


10.994 (P) 


10.972 


(P) 


18.995 


X6 2 (4P) 




10.963 


(P) 


11.020 (P) 


11.006 


(P) 


8.573 


Xb2(QF) 




10.931 


(P) 


11.196 (P) 


11.150 


(P) 


13.978 


Xb2($F) 




11.124 


(P) 


11.221 (P) 


11.187 


(P) 


18.995 


Xb2(QF) 




11.326 


(P) 


11.411 (P) 


11.380 


(P) 


13.978 , 



(83) 



and for 1 + states, hb(l + ): 



State 


Data I(bb) 


U{bb) 


R{bb) 


Rj 


h(lS) 


9.889 


9.902 


9.896 


0.922 


h b (2S) 


10.259 


10.255 


10.253 


2.782 


h b (3S) 


10.518 


10.530 


10.529 


5.781 


h b (AS) 


10.700 


10.768 


10.764 


18.839 


h(5S) 


10.759 


10.990 


10.971 


13.699 


h b (6S) 


10.921 


11.192 


11.200 


11.668 



(84) 



The wave functions for solution I(bb) are presented in Appendix II: we show the wave 
functions in the /c-representation and give coefficients q, f3 for Eq. (77). For 1 and 2 ++ 
states, we give Wll 1 '- one can see that the admixture of second components is small. 

3.2 Radiative decays (bb)i n —> 7(66) OM t 

Figure 3 shows the radiative transitions which are included into the fitting procedure. 



19 



The fit gives us the following values for radiative decay of T-mesons (partial widths in keV): 

Data 



Process 
T(1S) - 
T(2S) - 
T(2S) - 
T(2S) - 
T(2S) - 
T(2S) - 
T(3S) - 
T(3S) - 
T(3S) - 
T(3S) - 
T(3S) - 
T(3S) - 



^7] m {2S) 
!Xm{lP) 

!Xb2{lP) 

-frjbo(3S) 
7Xbo(2P) 
7X6i (2P) 
7X62 (2P) 



1.7±0.2 
3.0±0.5 
3.1±0.5 



1.4±0.2 
3.0±0.5 
3.0±0.5 



7(66) 


U{bb) 


R(bb) 


[7] 


[8] 


0.0092 


0.0100 


0.0079 






0.0025 


0.0015 


0.0008 


— 


— 


0.0006 


0.0002 


0.0002 


— 


— 


1.1023 


1.0669 


0.9611 


1.62 


1.41 


2.7288 


2.3675 


2.1268 


2.55 


2.27 


2.9100 


2.6674 


1.9367 


2.51 


2.24 


0.0016 


0.0007 


0.0005 






0.0013 


0.0000 


0.0002 






0.0006 


0.0001 


0.0001 






1.1191 


1.3746 


1.1253 


1.77 




3.4368 


4.0831 


3.1279 


2.88 




3.7266 


4.7438 


3.1686 


3.14 





(85) 



For the illustration, in (85) we present the results of Refs. [7, 8] . 

The radiative decay of XbJ is not included into fitting procedure. We have the following 
predictions for the partial widths (in keV): 



6 • 10- 2 
35 ±8) • 10~ 2 
22 ±4) • 10~ 2 
0.9 ±0.6) • 10~ 2 
4.6 ±2.1) • 10~ 2 
_,5± 1.3) • 10~ 2 
21.0 ±4.0) • 10~ 2 
7.1 ± 1.0) • 10~ 2 

16.2 ± 2.4) • 10- 2 16.14 18.25 14.41 (86) 



Process 




Data 


Xbo(lP) ~ 


+ 7 T(LS) 


< r 


tot(Xbo(lP) 


Xbi(lP) - 


- 7 T(15) 


r< i 


:xw(ip)) • 


Xm{lP) ~ 


-> 7 T(1S) 


r*ot 


:x6 2 (ip)) • 


Xm(2P) - 


+ 7 T(1S) 


r< i 


:x6o(2P)) • 


Xmi^P) ~ 


■+ 7 T(2S) 


r< i 


:x6o(2P)) • 


Xbi(2P) - 


- 7 T(15) 


r< i 


[Xbi(2P)) ■ 


Xbi(2P) - 


-> 7 T(2S) 




[Xbi(2P)) ■ 


Xb2(2P) - 


- 7 T(15) 




'^p)) ■ 


Xb2(2P) - 


-> 7 T(2S) 




:x b 2(2P)) ■ 


Xbo(3P) - 


+ 7 T(1S) 






Xbi^P) - 


+ 7 T(1S) 






Xb2(3P) - 


- 7 T(15) 






Xbo(3P) - 


-> 7 T(2S) 






XbifiP) - 


-> 7 T(2S) 






Xb2^P) ~ 


-> 7 T(2S) 






Xbo(3P) - 


-> 7 T(3S) 






Xbi(?P) - 


-> 7 T(3S) 






Xb2(3P) - 


-> 7 T(3S) 







I {lib) 


U(bb) 


R(bb) 


43.49 


52.79 


34.47 


53.31 


63.77 


42.62 


47.61 


56.15 


38.10 


7.49 


9.25 


7.08 


11.94 


15.88 


11.06 


12.76 


16.85 


11.69 


12.92 


14.40 


16.01 


18.11 


20.58 


12.37 


16.14 


18.25 


14.41 


2.69 


3.56 


2.84 


4.41 


7.06 


5.43 


5.60 


8.11 


6.01 


2.05 


1.86 


2.04 


3.38 


3.88 


4.09 


4.37 


4.59 


4.57 


7.59 


10.37 


7.19 


10.74 


15.85 


11.33 


9.90 


13.81 


9.87 



The total widths T tot (xbj(^P) and T tot (xbj(2P), with J = 0, 1,2, have not been measured yet. 

The calculations performed on the basis of (86) give us the following estimates for total 
widths: 



1 tot (Xbo(lP)) < 730 keV, 



(87) 
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r tot ( X 6i(lP)) ^ 120- 200 keV, 

r tot ( X62 (lP)) ^ 180- 270 keV, 

T tot ( X wC2P)) ~ 180-480 keV, 

r tot ( X 6i(2P)) ~ 50 - 80 keV, 

T tot (xb2(2P)) ^ 70 - 120 keV. 

The fit gives us the following values for partial widths of radiative decays of r/feo-mesons: 

Process Data 7(66) 

rj^S) -> 7 T(LS) 0.21 

77 60 (35) -> 7 T(1S) 0.27 

Vbo(3S) -> 7 T(2S) 0.05 



17(&6) 77(66) 

0.20 0.12 

0.18 0.11 

0.02 0.01 



(88) 



3.3 The 66 component of the photon wave function 

Fitting to the reactions with the 7 — > 66 transitions, we determined the parameters C n ,/3 7 , Sq 
defined in (63) for G S '° bl (s). For solutions 77(66), 7(66) and 17(66), they are as follows (in GeV): 



7(66) 




77(66) 


77(66) 


Cis ~- 


= -0.263 


Cis 


= -0.800 


Cis 


= -0.220 


C 2 s z 


= -0.296 


C2S 


= -0.303 


C2S 


= 0.092 


C3S z 


= 0.057 


C3S 


= 0.074 


C3S 


= -0.038 


Cis z 


= 0.298 


Cis 


= 0.197 


Cis 


= -0.027 




= -2.000 


C^s 


= -0.781 


Css 


= 0.262 


Ces ~- 


= 1.093 


Cqs 


= 2.000 


C6S 


= -0.441 


C\D ■ 


= -0.554 


C\d 


= -0.328 




= 0.168 


C2D : 


= -0.284 


C2D 


= 0.233 


C2D 


= 0.109 


6 7 = 


2.85 


6 7 = 


= 2.85 


6 7 = 


= 2.85 


s = 


18.79 


s = 


18.79 


s = 


18.79 



(89) 



Experimental values of partial widths included into fitting procedure as an input together with 
those obtained in the fitting procedure are shown below: 



Process 




Data 


7(66) 


U{bb) 


77(66) 


[39] 


T(LS) - 






1.314±0.029 


1.314 


1.313 


1.314 


1.01 


T(2S) - 






0.576±0.024 


0.576 


0.575 


0.576 


0.35 


T(3S) - 






0.476±0.076 


0.476 


0.476 


0.476 


0.25 


T(AS) - 






0.248±0.031 


0.248 


0.248 


0.248 


0.22 


T(5S) - 






0.31±0.07 


0.310 


0.310 


0.310 


0.18 


T(6S) - 






0.130±0.03 


0.130 


0.130 


0.130 


0.14 



The last column in (90) demonstrates the results of [39]. 

The vertices G^ bb -(s), which represent our solutions and given in (89), are shown shown 
in Fig 4. The data and predictions for the two-photon partial widths i] b0 — > 77, Xbo ~^ 77? 
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Xb2 77 are as follows: 
Process 

?7 6 o(l^) -> 77 
i] b0 (2S) -> 77 
776o(35') -> 77 
Xbo(l-P) -> 77 
Xbo(2P) -»• 77 
Xbo(3P) -> 77 
Xm(I-P) -> 77 
Xb2 (2P)^77 
X b2 (3P) -> 77 

(91) 

The last six columns represent the calculation results of [40-45]. 



Data 7(66) 


U(bb) 


R(bb) 


[40] 


[41] 


[42] 


[43] 


[44] 


[45] 


1.554 


1.851 


1.537 


0.35 


0.22 


0.46 


0.46 


0.45 


0.17 


1.928 


2.296 


1.906 


0.11 


— 


0.20 


0.21 


0.13 


— 


2.139 


2.547 


2.115 


0.10 


0.084 


— 


— 


— 


— 


0.024 


0.029 


0.021 


0.038 


0.024 


0.080 


0.043 






0.023 


0.028 


0.021 


0.029 


0.026 










0.023 


0.027 


0.020 














0.016 


0.020 


0.013 


0.0080 


0.0056 


0.0080 


0.0074 






0.015 


0.020 


0.013 


0.0060 


0.0068 










0.015 


0.019 


0.012 















3.4 Potentials in the solution I(bb) 

To be illustrative, let us demonstrate the interaction in one of the solutions, say, I(bb), in the 
language of potentials. 

In the solution I(bb), the confinement potential is due to scalar exchanges only, V^ f (r) = 
—0.151 + 0.160r (GeV), it is shown in Fig. 5a. For the illustration, we demonstrate also 
the ^levels created by this confinement potential alone. However, the scalar exchange has 
the short-range component v}^ rt (r) = 0.506 exp(— 0.2r) — 0.250/r/ea;p(— 0.2r) (GeV), see Fig. 
5b, which pushes the levels higher. The potential related to the t-channel vector exchanges, 
V( v \r) = 0.812-0.867exp(-0.4r)-0.300/r (GeV), does not contain the increasing part (~ r), 
see Fig. 5c. The last term may be interpreted as a one-gluon exchange, with a s = 0.300 ■ 3/4 ~ 
0.23. Also in Fig. 5c, we demonstrate for the illustration the ^levels, which would be created 
by vector exchange forces only. 

In the solution U(bb), the vector-exchange forces contain the one-gluon exchange term: 
Vshort( r ) = 1-355 exp(—0.5r) — 0.500/r (GeV) which corresponds to a s ~ 0.38. 



4 Conclusion 



In the framework of the method, which is in fact a variant of the dispersion relation approach, 
we have performed the description of the bottomonium spectra: the 66-levels and their radiative 
transitions such as (bb) in — > 7 + (bb) out , e + e~ — > V(bb) and bb-meson — > 77. Using quark- 
antiquark interaction as an input, we have obtained several variants of a reasonably good fit 
of the data. The ambiguities in a reconstruction of the soft region bb interaction underline the 
problem we face: a scarcity of the radiative decay data. To restore the bb interaction, one needs 
much more data, in particular, on the two-photon reactions: 77 — ■> bb-meson, including the 
bottomonium production by virtual photons in 77* and 7*7* collisions. 
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We pay a considerable attention to the presentation of the bb wave functions. The matter 
is that a solution may represent rather similarly the state levels, though with different wave 
functions. Therefore, we think that solution should be characterized by two characteristics, 
that is, the position of the level and its wave function. 
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5 Appendix I. The structure of pseudoscalar, scalar and 
vector exchanges 

The loop diagram, that includes the interaction, is given by the expression as follows: 

Sp[Q%f£(rn + fei)0/(m + K)Q^i\m - k)Oi(m - k 2 )] = v} s ' L ' J \-l) J 0^ , (92) 

where ki, k 2 are the momenta of particles before the interaction, k[, k' 2 are their momenta after 
the interaction, and the operators O/s are given by (14). 

For the singlet (S = 0) states in the case of scalar, pseudoscalar and vector exchanges, we 
obtain: 

Vp AJ) = ^(Azk - Am 2 - V^)k j Pj(z) , 
y(o,J,J) = ^(Azk + 4m 2 - \Gs>)k j 'Pj(z) , 

y{o,J,J) = ^(4^ -8m 2 )K J Pj(z) . (93) 

Here, Pj{z) are Legendre polynomials depending on the angle between final and initial particles 
and 

« = |k||k'|. (94) 

Near the threshold, the pseudoscalar interaction has a higher order factors k |k||k'| than in the 
case of scalar and vector interactions thus playing a minor role for mesons consisted of heavy 
quarks. The scalar and vector interactions in the lowest |k||k'| order are equal to each other in 
absolute value but have opposite sign. 

To obtain the expressions for triplet states, first, let us calculate the trace with vertex 
functions taken as 7^. Then, general expression can be obtained by the convolution of the trace 
operators: 

Spbv(ra + ki)Oi(m + fc' 1 )7 I/ (m - k' 2 )Oi(m - k 2 )] = (a[ + zKa{) g^ u +a{k^k^ + 
a 4^n k u + (a 5 + ZKa 6 ) k il k v + a 7 {k ll k u —k^ k u ) . (95) 
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The coefficients for the scalar, pseudoscalar and vector exchanges are equal to 



0/ 


1 


75 


7m 


a[ 


\ / ss'(Am 2 + \ / ss') 


\Jss'(Am 2 — yj ss') 


-2ss' 


a\ 


-A\fss' 


+4 v / ss 7 




a{ 


+As' 


-4s' 


-8s' 


a{ 


+As 


-4s 


-8s 


a{ 


4 (4m 2 - y/ss 7 ) 


4(4m 2 + v / ss 7 ) 


8(8m 2 — v/ss 7 ) 


a{ 


-16 


+16 


+32 


a\ 


+4 v / ss 7 


-Ay/ss 1 


+8 v / ss 7 • 



(96) 



Then, for S — 1 and L = J states we obtain: 



(i,j,j) 



ss'k j 



(Azn-Am 2 -\/s~s l )P J (z) 



An 



(zPj(z) - Pj^z)) 



JO 

V (1,J,J) = ^i K J 

7m 



J+l 

Azn-Am 2 + y/s~s 1 )P J {z) + -^-( z Pj(z) - Pj-i(z)) 



(2Vss' + 8zk)Pj{z) - 



J+l 



J+l 

(zPj(z) 



j -i 



») 



(97) 



Likewise, the states with L — J ± 1 are expressed as follows: 



(1,L,L',J) 



-1 



2J+ 1 



l+l' 
K 2 



/„,(£+') 



(98) 



fc=i 



We use additional index (L r ) to describe transitions between states with L + = J+l and L = J—l. 



(L,L>) 

(L,L r ) 
V 2 

V (L,L>) 

JL,L>) 
V 4 

JL,L>) 
v 5 

JL,V) 
u 6 
{L,L') 



(2J+l)P J+1 (z) 

(2J+1)zkP j+1 (z) 

-(J+l)P J+1 (z)\k\ 2 

-(J+l)P J+1 (z)\k>\ 2 

-{J+1)kPj{z) 

-(J+1)zk 2 Pj(z) 

(2J+1)k(zP j+1 (z)-Pj(z)) 



■L+ 



L+- 



(2J+l)Pj^(z) 
(2J+1)zkP J - 1 (z) 
-JPj^(z)\k\ 2 
-J Pj^(z)\k'\ 2 
-JkPj(z) 

-JzK 2 Pj(z) 






AkPj+i 

APj(z)\k'\ 2 

AzKPj(z)\k'\ 2 





,JeI1 

K 



AkP j+1 (z) 
APj(z)\k\ 2 
AzKPj(z)\k\ 
. 



(99) 



Here, A = \JJ(J + 1) and k are defined by (94). 

6 Appendix II. Wave functions in the bb sector 

The tables 1-4 demonstrate the coefficient values c\ n \ which determine the wave functions 
^j(s,l,j) accorc lj n g to (77). In the tables we also show Wlu, which determine the normalization 
condition, see (13). In Figs. 6-9, we demonstrate these wave functions. 
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Table 1: Constants c) n) from (77) (in GeV) for the wave functions of the T-mesons in solution 
I(bb); also we present the normalization coefficients W 00: W 02: W 22 given by (12), (13). 





T(1S) 


T(2S) 


T(3S) 




Woo 




W 02 


w 22 


Woo 




w 02 


w 22 


Woo W 02 


w 22 




1.00050 - 


0.00060 


0.00010 


1.00101 - 


0.00122 


0.00021 


1.00115 -0.00143 


0.00028 


i 


^(1,0,1) 


^,(1,2,1) 


^(1,0,1) 


^(1,2,1) 


^,(1,0,1) 


^,(1,2,1) 


1 


2.0060 






-1.9545 


2.7111 






1.7177 


-4.7273 




1.9830 


2 


-10.8380 






13.5543 


9.8323 






-12.2344 


9.5837 




-12.8048 


3 


73.0097 






-36.0805 


-72.4706 






32.3653 


-52.0435 




33.3246 


4 


-230.4699 






47.0494 


221.2708 






-41.5308 


240.2794 




-42.7176 


5 


397.0086 






-30.8840 


-386.7066 






26.4284 


-457.6044 




26.8052 


6 


-389.9981 






7.9660 


386.9034 






-6.1016 


453.3999 




-5.5935 


7 


218.1722 






1.2842 


-219.7120 






-1.5700 


-251.5948 




-2.0591 


8 


-64.5122 






-1.1142 


65.7237 






1.0735 


73.9804 




1.2288 


9 


7.8476 






0.1657 


-8.0769 






-0.1519 


-9.0153 




-0.1688 




T(4S) 


T(5S) 


T(6S) 




W 00 






w 22 






w 02 


w 22 


Woo W 02 


W 22 




1.00118 - 


0.00155 


0.00037 


1.00149 - 


0.00185 


0.00036 


0.80131 0.03464 


0.16405 


i 




^,■2,1) 


^(1,0,1) 






^,(1,2,1) 


1 


-4.6053 






-1.9248 


4.7469 






-2.3093 


4.3499 




-33.8613 


2 


-12.9581 






16.3285 


26.3926 






11.4831 


66.6334 




71.5349 


3 


122.1101 






-45.3030 


-218.0913 






-23.7269 


-672.1165 




95.8832 


4 


-232.1932 






57.4674 


306.9638 






26.0199 


2023.7316 




-345.3417 


5 


242.2475 






-34.4863 


140.7892 






-15.4406 


-2851.3726 




273.7456 


6 


-211.8979 






6.2974 


-597.4634 






3.9824 


2118.8719 




-24.6291 


7 


139.6440 






2.9879 


479.8989 






0.3477 


-850.6343 




-61.8950 


8 


-50.5914 






-1.5944 


-161.4547 






-0.4002 


175.5466 




27.9503 


9 


7.2976 






0.2100 


20.2973 






0.0570 


-15.0191 




-3.6706 
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Table 2: Constants c\ n> from eq. (77) (in GeV units) for the wave functions of T-mesons in the 
solution I(bb); normalization coefficients W 00: W 02: W 2 % are given by (12), (13). 





T(1D) 


T(2L>) 


T(3L>) 




Woo 


W 02 


w 22 


Vt^oo 




Wo 2 


w 22 


Woo W 02 


w 22 




0.00083 -0.00214 


1.00131 


0.00144 - 


0.00360 1.00217 


0.00210 -0.00499 


1.00290 


i 




^,(1,0,1) 


^(1,2,1) 


^(1,0,1) 




^,(1,0,1) 


1 


0.5316 




0.0051 


-14.2711 






-0.0804 


6.0089 




0.0734 


2 


22.2640 




0.2961 


29.5676 






0.4676 


104.3726 




0.4895 


3 


-61.8266 




-2.1943 


-62.5673 






-3.1827 


-367.7237 




-4.4838 


4 


79.4272 




7.1306 


100.0037 






11.9276 


446.4882 




12.5128 


5 


-50.4830 




-12.7875 


-73.1973 






-21.9636 


-236.8819 




-20.8090 


6 


10.6083 




12.9455 


15.1171 






22.1691 


31.8558 




22.5104 


7 


4.5874 




-7.4078 


8.7956 






-12.6656 


23.2527 




-14.3861 


8 


-2.8530 




2.2290 


-5.1628 






3.8222 


-10.8663 




4.7718 


9 


0.4300 




-0.2751 


0.7691 






-0.4749 


1.4288 




-0.6355 




T(4L>) 


T(5L>) 


T(6D) 




W 00 


Wo 2 


w 22 


Woo 




w 02 


w 22 


Woo W 02 


W 22 




0.00230 -0.00585 


1.00355 


0.20433 - 


0.04594 0.84162 


0.00126 -0.00134 


1.00008 


i 


^(1,2,1) 








^,(1,2,1) 




1 


-32.4189 




-0.0913 


80.9176 






2.0418 


-377.8732 




0.5165 


2 


-32.2563 




-0.3537 


-211.1041 






33.2603 


2121.7911 




-1.1281 


3 


385.0197 




1.9186 


-54.4316 






-333.9307 


-4529.4955 




-14.6612 


4 


-592.6611 




6.0232 


550.7145 






1012.1198 


4646.2765 




69.2347 


5 


321.7644 




-30.8630 


-476.3406 






1440.4555 


-2245.1955 




123.0232 


6 


-16.3314 




44.9067 


41.8114 






1084.8218 


251.1282 




109.8199 


7 


-45.1952 




-29.9473 


113.3274 






-443.3884 


204.4454 




-52.1311 


8 


15.1386 




9.5444 


-50.9272 






93.7641 


-81.9575 




12.5435 


9 


-1.4251 




-1.1922 


6.6630 






-8.2835 


9.2592 




-1.2099 
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Table 3: Constants c\ n> from (77) (in GeV) for the wave functions of Xfc2-mesons in the solution 
I(bb) normalization coefficients W q,W 02 ,W22 are given by (12), (13). 





X62(1P) 


Xb2(2P) 


Xb2(3P) 




Wii ^13 


w 33 


Wn W 13 


w 33 


Wn W 13 


w 33 




1.00075 -0.00094 


0.00019 


1.00128 -0.00203 


0.00075 


0.98981 -0.00494 


0.01514 


i 


^(1,1,2) 


^,(1,3,2) 


</,(l,l,2) 


^,(1,3,2) 


^,(1,1,2) 


^,(1,3,2) 


1 


13.0540 




-4.1772 


-34.6278 




16.0313 


-257.7504 




100.2644 


2 


-114.9374 




23.4967 


456.0158 




-90.5844 


2844.0283 




561.9424 


3 


487.8876 




-52.5492 


-1939.7651 




201.0440 


-12019.5471 


1240.5753 


4 


-1062.2386 




60.1441 


4168.0416 




227.9811 


25892.6834 


-1399.5333 


5 


1317.9252 




-37.5852 


-5113.9135 




141.4654 


-31653.5893 




863.3132 


6 


-967.3832 




12.2496 


3712.1978 




-46.2459 


22808.4240 




280.4788 


7 


414.8491 




-1.5240 


-1571.7829 




6.1338 


-9564.2325 




37.0697 


8 


-95.8539 




-0.1252 


357.7836 




0.2835 


2152.2305 




1.6274 


9 


9.2113 




0.0356 


-33.7704 




-0.1074 


-200.3714 




-0.6277 




X62(4P) 


X62(5P) 


Xb2(6P) 




W n W 13 


W 33 


Wn W 13 


w 33 


Wn W 13 


w 33 




1.00071 -0.00246 


0.00176 


1.00246 -0.00318 


0.00072 


1.00311 -0.00397 


0.00086 


i 


^1,1,2; 


^1,3,2) 


^,(1,1,2) 


^,(1,3,2, 


^(M,2) 


^,(1,3,2, 


1 


61.1310 




-27.9538 


-21.7721 




4.7026 


10.8365 




-2.3155 


2 


-804.3179 




157.6290 


68.2854 




-23.9703 


144.6924 




5.6067 


3 


3433.8802 




347.3610 


-194.8198 




48.5937 


-929.0568 




2.8923 


4 


-7271.1541 




388.8039 


630.5268 




-51.7634 


1868.8590 




-20.3080 


5 


8755.3125 




236.5028 


-1080.4598 




32.0724 


-1686.1180 




23.4272 


6 


-6242.9707 




74.8939 


933.0968 




-12.0176 


679.5364 




-11.5459 


7 


2595.9632 




-9.1414 


-421.0616 




2.7404 


-61.8407 




2.2550 


8 


-579.2785 




-0.6539 


94.9352 




-0.3698 


-31.7481 




0.0334 


9 


53.3862 




0.1874 


-8.3610 




0.0242 


6.7593 




-0.0471 




Xb2(lF) 


Xb2(2F) 


X62(3F) 




Wn W 13 


w 33 


Wn W 13 


w 33 


Wn W 13 


w 33 




0.00058 -0.00191 


1.00132 


0.00097 -0.00315 


1.00218 


0.00124 -0.00390 


1.00266 


i 


0(1,3,2, 


^,(1,1,2) 


0(1,3,2, 


^,(1,1,2) 


^,(1,3,2) 


^,(1,1,2) 


1 


-3.6248 




-0.0963 


-14.7982 




0.0560 


-37.1706 




2.9107 


2 


-7.2244 




0.8055 


0.1033 




-1.0323 


9.5966 




-34.8030 


3 


19.9292 




-2.9097 


31.5118 




3.4274 


147.9125 




151.9065 


4 


-17.8259 




5.3091 


-19.0518 




-4.1874 


-177.1377 




332.1606 


5 


5.4336 




-5.0342 


-2.9800 




1.6200 


36.1365 




410.3068 


6 


2.6649 




2.4892 


3.3342 




1.1027 


48.9104 




298.9169 


7 


-2.9157 




-0.5206 


1.0840 




-1.4679 


-35.1616 




126.9788 


8 


0.9614 




-0.0152 


-0.9814 




0.5902 


9.1555 




-29.0017 


9 


-0.1144 




0.0160 


0.1614 




-0.0846 


-0.8818 




2.7474 




Xb2^F) 


Xb2(5F) 


Xb2(6F) 




W n W 13 


w 33 


Wn W 13 


w 33 


Wn W 13 


w 33 




0.00244 -0.00497 


1.00253 


0.00119 -0.00465 


1.00346 


0.20223 0.04827 


0.74950 


i 


0(1,3,2) 


W,(l.l,2) 


^,(1,3,2) 


^,(1.1.2) 


W,(l>3,2) 


.0(M,2) 


1 


137.2096 




-3.9590 


-295.5265 




0.0318 


-342.8625 




88.3149 


2 


-395.5970 




55.7513 


1172.2960 




-8.1102 


1605.5381 




670.5111 


3 


346.9297 




259.4898 


-1660.4991 




60.4182 


-2779.6979 


1944.5425 


4 


-64.2683 




578.6142 


946.4669 




176.4908 


2084.3230 


-2837.7766 


5 


-14.0966 




709.0203 


-65.5667 




263.0556 


-374.7669 


2282.8464 


6 


-40.3793 




502.8571 


-139.7540 




217.9552 


-406.6671 


-1033.5668 


7 


38.2866 




205.4367 


44.294^ 


7 


101.2739 


272.3230 




257.4060 


8 


-11.2262 




44.7930 


-0.5172T 




-24.6525 


-64.5095 




-33.5769 


9 


1.0954 




-4.0251 


-0.8887 




2.4484 


5.5083 




2.1172 



Table 4: Constants c\ from (85) (in GeV) for wave functions rjbo, Xbo, Xbi an d , i n solution 
I(bb) 



i 


IU\) \ / 


Vbo(2P) 


Vbo(3P) 


%o(4P) 

luvj \ / 


Vbo(5P) 


Vbo(QP) 


1 


4.5338 


-1.1326 


11.9074 


-18.3892 


29.5281 


-3.5400 


2 


-55.4035 


73.3964 


-136.8130 


411.7847 


-480.8391 


201.4420 


3 


329.9544 


-440.1212 


798.7926 


-2495.6136 


2863.6920 


-1321.6947 


4 


-942.3382 


1240.9035 


-2311.1680 


6933.5200 


-8155.6641 


3472.8097 


5 


1474.8959 


-1923.5277 


3573.4205 


-10437.4060 


12468.2246 


-4685.6254 


6 


-1326.7453 


1712.0637 


-3132.4041 


9006.4575 


-10761.0759 


3555.8036 


7 


682.8131 


-870.3125 


1560.9990 


-4427.9932 


5234.8626 


-1537.4724 


8 


-186.4237 


234.3339 


-411.2328 


1148.6087 


-1334.8254 


352.9427 


9 


20.9478 


-25.9307 


44.4414 


-121.6868 


138.3033 


-33.1141 


i 


XboCi-P) 

/\.U\J \ / 


Xbo(2P) 


/V I'U \ J 


/V I'U \ / 


A, u\) \ / 


Xw(6P) 


1 


29.2420 


-61.6937 


-219.3775 


126.2557 


-42.4376 


-6.8271 


2 


-304.0472 


755.4596 


2411.7486 


-1536.0712 


316.4890 


347.2500 


3 


1304.2479 


-3236.5223 


-10285.1211 


6586.8209 


-1275.2749 


-1772.0896 


4 


-2863.9307 


7047.5098 


22382.0814 


-14186.9213 


2956.9372 


3636.4974 


5 


3585.8612 


-8744.9017 


-27622.3384 


17341.7384 


-3876.4948 


-3793.2353 


6 


-2656.7748 


6413.5855 


20084.3830 


-12510.0465 


2897.0308 


2166.3382 


7 


1149.8200 


-2743.1742 


-8497.3906 


5252.0887 


-1220.4556 


-673.3240 


8 


-268.0686 


630.9120 


1929.4613 


-1182.3208 


268.9864 


102.7608 


9 


25.9818 


-60.1922 


-181.3145 


109.9412 


-23.9653 


-5.3795 


i 






Xbi&P) 


Xbii^P) 


Xbi(5P) 


Xbi(QP) 


1 


33.0397 


239.7994 


-164.4051 


34.1495 


26.9515 


8.4967 


2 


-332.6519 


-2634.9469 


1721.8106 


-485.7810 


-126.6739 


149.8510 


3 


1385.1827 


10912.2385 


-7107.3970 


2029.1900 


430.9352 


-885.8464 


4 


-2958.4499 


-23096.8480 


15028.0903 


-4151.1836 


-1094.3863 


1678.5809 


5 


3605.9605 


27851.8047 


-18066.4365 


4843.3153 


1580.6111 


-1375.8433 


6 


-2602.3368 


-19853.0333 


12821.3870 


-3359.3446 


-1243.8826 


421.1791 


7 


1097.2744 


8252.3736 


-5300.9385 


1361.1611 


531.5536 


55.6201 


8 


-249.2044 


-1843.7114 


1176.8983 


-295.8920 


-115.6049 


-59.6685 


9 


23.5136 


170.6412 


-108.1134 


26.4999 


9.9105 


9.4863 


i 


h b (lP) 


h b (2P) 


h b (3P) 


h b {AP) 


h b (5P) 


h b (6P) 


1 


-26.9406 


245.4519 


-158.2896 


29.3812 


24.3700 


10.9592 


2 


265.1898 


-2696.0191 


1654.8622 


-435.3349 


-96.8333 


123.8420 


3 


-1104.2180 


11165.1311 


-6832.1653 


1822.9779 


306.2952 


-785.0079 


4 


2358.6741 


-23627.7741 


14450.3845 


-3718.3930 


-838.9663 


1479.3562 


5 


-2874.8697 


28482.8717 


-17375.2937 


4327.6261 


1290.5769 


-1149.1916 


6 


2074.5354 


-20294.0610 


12332.5098 


-2998.2763 


-1052.4158 


266.7467 


7 


-874.6420 


8431.1444 


-5099.3191 


1214.5985 


458.5851 


117.4293 


8 


198.6257 


-1882.4288 


1132.2297 


-264.0975 


-100.7958 


-72.9291 


9 


-18.7411 


174.0887 


-104.0177 


23.6613 


8.6824 


10.6504 
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(0 (g) (h) 

Figure 1: a,b) Diagrams for the two-photon decay of QQ state and c) cuttings in the spectral 
integral representation. Initial (d) and final (e) state interactions of quarks in the decay dia- 
grams, f) Graphical representation of the spectral integral equation for the QQ vertex. g,h) 
Interaction of quarks in the vertex QQ — > 7 and its approximation by the sum of transitions 
QQ - E V(n) - 7 . 

n 
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+ 




Figure 2: Quark transition diagram for the process e + e — > V . 




Figure 3: Radiative decays of the bottomonium systems, which were taken into account in 
the fit (solid lines). The dashed lines show radiative transitions with the known ratios for the 
branchings Br[xbj{^P) — > 7T(2S , )]/.Br[x&j(2P) — > 7T(1S')] , these ratios are not included in 
the fit. 
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1.5 2 2.5 3 3.5 4 4.5 
k, GeV 




Figure 4: a) Vertices G^_^ (a) and G^J^ (b) for the solutions I(bb) (solid line), U(bb) (dash 
line) and R(bb) (dotted line). 
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l(bb)-solution 



> ° 

a> 

O 

S 2 
> 

1 

-1 
-2 
-3 



V( s o )(r)=-0.151 + 0.160r 



J i i i I i i i I i i i I i i i I i i i L 



8 



10 12 
r GeV 1 



> ° 

a> 

O 

S 2 
> 

1 

-1 




- 

- 
- 


b) 








V( s h ) ort (r)=0.506e 02r -°- 2 r 50 e 02r 


■ i i I i 


i 1 i i i 1 i i i 1 i i i 1 i i i 1 i 


) 2 


4 6 8 10 12 




r GeV" 1 



m O I i I i I i i i I i i i I i i i I i i i I i i i I i i i I i i i I i i i I i i 

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

r GeV" 1 

Figure 5: Potential for the solution I(bb): confinement potential (a), short-range repulsive one 
(b) and attractive potential due to vector exchange (c). To be illustrative, we show the CT + 
levels, which would be created by the attractive potentials. 
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Figure 6: Wave functions for T(nS) in the solution I(bb). Solid and dashed lines stand for 
^,(i,o,i) an d ^,(i,2,i) . 
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Figure 7: Wave functions for T(nD) in the solution I(bb). Solid and dashed lines stand for 

^(1,0,1) and ^,(1,2,1). 
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Figure 9: Wave functions for Xbo in the solution 7(66). 
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Figure 11: Wave functions for Xb2 in the solution I(bb). Solid and dashed lines stand for ■0( 1 ' 1 ' 2 ) 
and^ (1 ' 3 ' 2) . 
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Figure 12: Wave functions for Xb2 in Solution 7(66). Solid and dashed lines stand for tp^ 1 ' 1 
and V (1,1 ' 2) - 



42 




43 



